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ABSTRACT 


Let S: pxp be distributed as central Wishart matrix with n degrees of freedom 
and let E(S) = lb where I, is an identity matrix. Also, let 0) <- +» <6,, be the 
roots of S. In this report, the authors gave exact upper 10%, 5%, 2.5% and 1% 
points of the distributions of 6; (i=1,2,...,p-1) for p=2(1)10 and n=(p+1) (1) 20 


(2) 30(5)50. 


iii 


TABLE OF CONTENTS 


Section Page 
1. Introduction 1 
2. Cumulative Distribution Yunction of the Individual 1 
Roots 
53. References 5 


Tables 


iv 


1. INTRODUCTION 


The marginal distribution of the individual roots of the Wishart matrix 
are useful in testing certain statistical hypotheses, Pillai and Chang [7] 
constructed tables for the upper percentage points of the largest root of the 
Wishart matrix. Davis [2] gave a recurrence relation for the marginal densities 
of the individual roots by using results in [1]. Krishnaiah and Waikar [5] gave 
exact expressions for the cumulative distribution functions (c.d.f.'s) of inter- 
mediate roots. In this paper, we give exact upper 10%, 5%, 2.5% and 1% points 


of the smallest and intermediate roots of the Wishart matrix. 


2, CUMULATIVE DISTRIBUTION FUNCTIONS OF THE INDIVIDUAL ROOTS 


Let S:pxp be distributed as Wishart matrix with n degrees of freedom. and 
let E(S) =n Ty? where 1, is the pth order identity matrix. Also, let Oy<eees ey 
be the eigenvalues of S. Then the joint density of 04 <++- <6, is given by 


£ = k n fet 4 
(01 4+++ 58.) z Oe exp (- as )] 7 (85-65) 


1>j 
0<64<+ +640, <0 (2.1) 
where 
r= (n-p-l) /2, and 
n 
P 
k(pin) = w/? (4) EF eran /yrC-0 7291, (2.1a). 
1= 


The following exact expression for the'c.d.f. of the intermediate root 6, (1<s<p-1) 
was given in Krishnaiah and Waikar [5]: 
Plo.<x] = Ple,,4<x] + k(p.n)},#0(w35,{k,,.-..k,}, 05x). 


e(p; pes, {tyseees t X,~), (2.2) 


pes? 


where {kyy.++,kg) is a subset of integers {0,1,...,p-1} such that ky<.e eck, 

and ty<s+.cth is the subsct complementary to k,<...<k, while }, denotes the 
summation over (°) possible subsets ky<ee esky. Further w{y) = exp (- » ) and the 
sign inside by is positive or negative according as s(s+3)/2 + 5 K, is even or odd. 


i=] 
The function p(-+) is defined by 


pv; Py {ky y+. 5k), L,U) = ACY; am, {ky ,+++ sk) LU) when p = 2m (2.5) 


and 
NEN ote 
e(Ws Detkyse sk ds TU) = 2 CL Ty LU). 
i=l i 
Gi(¥5 21, ky, 2.6 kong )s IeyU) when p = 2m+l (2.4) 


where 
3 


kK, wt 1/2 
4(y; am, {ky,---rko id, L,U) = CF, G4) i,j=l, yam], 


CLC) 2m Ky soang Key) dsU) 


“(6 (I.t)) i,jet,..., tel, tel,...,ame1|?/? 
1 
for t=1,...,2m+1 while Gi 1,k,,L,U)=1. Further 
t ~ rt rT 
fe (LU)= FS (L,U)-FE (LU), 
t U t 
re(L,u) = § F_ (1,8) e° w(0) de, 
S$ t Ss 


9) 
r. (L,8) -¢ x5 y(x) dx 
Also it is known (see[3,4]),that 
Ple,sx] = 1-k(p,n)e(¥sdir,c,=), (25) 


Jsing (2.2) and (2.5), we have constructed tables for the exact values of x for 


922(1)10, n=(p+1) (1) 20(2) 30(5)50,k=1(1) (p-1) and %0.01, 0.025, 0.05, 0.10 where 


Plo. é x) © l-g, 


The values are given up to four decimal places and they may differ from actua] 


values by at most one «mit in the last decimal. As a check for the accuracy of 


the tables we have used the programs to compute the values of 


Pixso5] + P[6, <xso J+... 


$6, + PLO, <xs8,] + P[0 <x] 


for some x and found them to differ 


fron 1 in the 12th decimal only, 


[6] 


[7} 
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The entries in the. following table are the values of x for different value 
‘ \ 

of N,?,S and « where 

Plagex] = (1-a), 


Dp 
aoe Pete . 
F(O4.+++ Op) = ee (0. il 0,/2)] a 
\ ; ‘ 


N<ay<. ; eae os SAycm , 


r= (N-P-1)/2, and k(P,\N) is given by fq. (2.14). ‘ 
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